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PREFACE 


The  Seismic  Hazards  Mapping  Act  (Chapter  7.8,  Sections 
2690  et  seq.,  California  Public  Resources  Code)    requires 
the  State  Geologist,  chief  of  the  Department  of 
Conservation's  Division  of  Mines  and  Geology  (DMG),  to 
designate  seismic  hazard  zones.  These  zones  assist  cities 
and  counties  in  fulfilling  their  responsibilities  for  protecting 
the  public  health  and  safety  from  the  effects  of  strong 
ground  shaking,  earthquake-induced  landslides,  liquefac- 
tion, or  other  ground  failures.  To  assist  the  State  Geologist 
in  fulfilling  this  responsibility,  the  Act  directs  the  State  Min- 
ing and  Geology  Board  (SMGB),  in  consultation  with  an 
advisory  board,  to  develop  guidelines  and  criteria  for  the 
preparation  of  seismic  hazard  zones  in  the  state.  This  re- 
port presents  the  recommendations  of  the  Seismic  Hazard 
Mapping  Act  Advisory  Committee  as  accepted  by  the 
SMGB.  It  is  expected  these  criteria  will  continue  to  evolve 
as  our  understanding  of  seismic  phenomena  and  the 
methods  used  to  assess  their  likelihood  and  potential  im- 
pacts on  the  built  environment  improve. 

The  Seismic  Hazard  Mapping  Act  Advisory  Committee 
formed  three  working  groups  composed  of  acknowledged 
experts  to  address  ground  shaking,  liquefaction,  and  land- 
slide hazards  in  an  effort  to  gain  a  consensus  on  how  to 
prepare  the  various  maps  (see  Acknowledgments). 


A  fourth  working  group  on  planning  and  implementa- 
tion was  formed  to  ensure  that  the  resulting  seismic 
zone  maps  would  be  of  practical  use  in  the  local  plan- 
ning and  building  department  decision-making.  Rec- 
ommendations from  these  working  groups  are  principal 
components  of  this  document. 

The  previous  unpublished  version  of  this  publication 
presented  criteria  for  delineating  liquefaction  hazard 
zones  and  recommended  that  current  methods  of 
evaluating  earthquake-induced  landslides  be  investi- 
gated. Furthermore,  because  the  potential  for  amplified 
ground  shaking  cannot  be  estimated  with  sufficient  reli- 
ability, the  previous  version  recommended  that  such 
hazard  zones  not  be  established  until  the  new  defini- 
tions for  site  conditions  are  released  in  the  Uniform 
Building  Code  (UBC)  (ICBO,  1997).  A  basis  for  delin- 
eating earthquake-induced  landslides  has  now  been 
adopted,  and  constitutes  the  principal  change  in  this 
document.  Although  new  definitions  of  site  factors 
have  been  adopted  in  the  UBC,  a  consensus  has  yet  to 
develop  by  the  Seismic  Hazard  Mapping  Act  Advisory 
Committee  on  whether  to  establish  hazard  zones  for 
amplified  shaking.  Decisions  in  that  regard  may  form 
the  basis  for  an  update  of  this  document. 


VI 
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PROBABILISTIC  SEISMIC  HAZARD  MAP 


INTRODUCTION 

The  California  Department  of  Conservation's  Division  of 
Mines  and  Geology  (DMG)  is  charged  with  implementing 
requirements  of  the  Seismic  Hazards  Mapping  Act  of 
1990.  Appropriate  maps  of  expected  ground  shaking 
hazard  are  required  and  are  an  underpinning  for  map- 
ping any  and  all  seismic  hazard  zones.  The  following 
recommendations  are  provided  to  assist  DMG  in  map- 
ping ground  shaking  hazard  on  a  regional  scale 
throughout  the  state. 

General  Considerations  for  Mapping 
Expected  Ground  Shaking  Hazard 

The  Advisory  Committee  recommends  preparation  of  a 
suite  of  regional  ground  shaking  hazard  maps  using 
Probabilistic  Seismic  Hazard  Analysis  (PSHA)  techniques 
(NRC,  1988).  The  following  maps  should  be  produced  at 
statewide  scales: 

1 .  Maps  of  peak  ground  acceleration  and  spectral  accel- 
eration at  0.3  sec,  1.0  sec,  and  3.0  sec,  with  exceedance 
probabilities  of  10%  in  50  years,  50%  in  50  years,  and 
10%  in  100  years. 

2.  Maps  of  peak  ground  acceleration,  weighted  with  re- 
spect to  a  M7.5  earthquake,  for  evaluation  of  liquefaction 
potential  and  earthquake-induced  landslide  potential, 
with  exceedance  probabilities  of  10%  in  50  years,  50%  in 
50  years,  and  10%  in  100  years. 

Existing  probabilistic  seismic  hazard  computational  codes 
are  acceptable  and  no  basic  modeling  developments  nor 
substantive  computational  code  changes  are  needed. 
The  results  should  capture  and  display  uncertainties 
about  input  parameters,  including  seismic  sources,  earth- 
quake frequency,  maximum  magnitude,  seismic  wave 
attenuation,  and  site  response.  Input  parameters  should 
be  developed  by  consensus  of  an  earth  science  team  us- 
ing consistent  approaches  throughout  the  state  and 
formal  uncertainty  elicitation  procedures  (NRC,  1977). 

PSHA  mapping  should  extend  to  the  near  offshore  re- 
gions and  use  Uniform  Building  Code  soft  rock 
conditions  as  the  base  site  condition  and  reference  soil 


column.  A  companion  report  should  be  prepared  that 
analyzes  the  key  sources  of  uncertainty  in  enough  depth 
and  detail  to  permit  users  to  factor  uncertainty  into  their 
use  of  the  maps.  The  analysis  of  uncertainty  may  require 
modest  computational  code  development.  Work  should 
be  coordinated  with  ongoing  PSHA  efforts  of  the  U.S. 
Geological  Survey  (USGS). 

Seismic  Source  Modeling 

Three  general  types  of  seismic  sources  are  expected, 
1)  sources  that  model  active  faults,  2)  sources  that  model 
"active"  structures  that  may  contain  significant  faults 
(i.e.,  active  fold  belts,  such  as  those  along  the  western 
edge  of  the  Central  Valley  and  within  the  LA  Basin),  and 
3)  sources  that  model  distributed  seismicity  that  cannot 
be  assigned  to  specific  geologic  structures.  All  three 
types  of  sources  can  be  readily  modeled  within  existing 
computational  programs.  The  details  of  fault  geometry 
should  not  have  a  major  impact  on  the  results  of  a  re- 
gional hazard  study  in  terms  of  its  effect  on  the  density 
function  for  distance  to  rupture.  (It  may  have  a  signifi- 
cant impact  on  parameters  such  as  maximum  magnitude 
and  seismicity  rate,  if  moment  (slip)  rate  methods  are 
used.)  Some  special  attention  to  details  of  geometry 
may  be  needed  in  the  northwest  to  model  the  Cascadia 
subduction  zone. 

The  seismic  sources  can  be  identified  on  the  basis  of  ex- 
isting extensive  fault  mapping  and  surface  and/or 
subsurface  mapping  of  actively  deforming  folds  for  Cali- 
fornia. Careful  thought  needs  to  be  given  to  "back- 
ground" sources  to  account  for  possible  unidentified  ma- 
jor sources.  Uncertainty  in  sources  can  be  modeled  by 
providing  weighted  alternatives. 

Maximum  Magnitude 

Maximum  magnitudes  for  fault-specific  sources  should 
be  based  on  interpretations  of  the  potential  maximum 
size  of  rupture  and  the  well-developed  empirical  rela- 
tionships between  rupture  dimensions  and  magnitude 
that  are  documented  in  the  literature.  Assessments  of 
maximum  magnitudes  for  tectonic  structures  may  have 
to  rely  more  on  analogy  than  on  specific  dimensions  of 
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structures,  although  the  general  characteristics  of  the 
structure  (e.g.,  long  and  continuous  folds  versus  short 
and  offset  folds)  may  suggest  trends  in  the  maximum 
size  that  could  be  used  to  weight  the  various  analogies. 
Assessments  for  seismicity  zones  and  background  zones 
most  likely  will  have  to  rely  on  arguments  based  on 
analogy,  largest  observed  events  without  surface  rupture 
manifestations,  and  historical  observations.  Uncertainty 
on  maximum  magnitude  should  be  modeled  using  a 
variable  with  a  distribution  rather  than  a  single  value. 

Earthquake  Frequency 

The  primary  model  for  earthquake  recurrence  should  be 
the  Poisson  model,  because  we  know  little  more  than 
average  rates  for  the  vast  majority  of  seismic  sources. 
Time-dependent  models  may  be  applicable  in  a  few 
areas.  This  could  be  tested  to  assess  how  regional  map- 
ping results  might  be  adjusted.  For  fault-specific  sourc- 
es, earthquake  frequency  (slip  rate)  should  be  based 
primarily  on  geologic  information  for  those  faults  where 
data  on  paleoseismicity  can  be  used  to  establish  a  rate. 
For  other  tectonic  structures,  other  geologic  information 
may  have  some  use  in  areas  where  rates  of  deformation 
can  be  established  and  where  a  fraction  can  be  attribut- 
ed to  movement  on  faults.  However,  historical 
seismicity  rates  will  likely  be  the  primary  source  of  recur- 
rence information  for  these  other  structures,  as  it  will  be 
for  distributed  seismicity  zones.  Recurrence  parameters 
should  be  modeled  as  variables  with  distributions. 


Minimum  Magnitude 

It  is  recommended  that  the  minimum  magnitude  of  in- 
terest be  set  about  M5.  It  may  be  desirable  to  compute 
results  for  a  higher  minimum  magnitude  to  capture  the 
level  of  hazard  from  major  earthquakes  compared  to 
the  hazard  from  moderate  earthquakes. 

Seismic  Wave  Attenuation 

A  new  generation  of  seismic  wave  attenuation  curves 
should  be  developed  using  an  updated  empirical  data- 
base from  recent  strong-motion  recordings.  This  work 
should  be  coordinated  with  ongoing  seismic  wave  atten- 
uation studies  at  the  USGS.  "Standard"  attenuation 
curves  should  be  developed  for  various  UBC  site  soil 
conditions. 

Magnitude  dependence  of  attenuation  dispersion 
should  be  confirmed  and  incorporated  into  the  PSHA  if 
appropriate. 

A  number  of  site/source/path  conditions  may  influence 
seismic  wave  attenuation.  Not  all  of  these  conditions 
are  accommodated  in  the  empirical  curves  when  they 
are  applied  at  a  given  site  (e.g.,  long  period  ground  mo- 
tions in  basins,  faulting  style,  near-source  effects  at  long 
periods,  crustal  structure,  focal  depth  and  topography). 
The  PSHA  should  proceed  with  an  awareness  of  these 
effects  and  they  should  be  discussed  in  the  commen- 
tary. In  general,  until  more  definitive  procedures  can  be 
developed,  the  PSHA  should  treat  these  effects  as  part 
of  the  randomness  in  seismic  wave  attenuation. 


AMPLIFIED  SHAKING  HAZARD  ZONES 


Building  codes  are  currently  the  primary  means  of  miti- 
gating the  effects  of  strong  earthquake  shaking  on 
buildings.  The  effect  of  local  surface  geology  on  expected 
shaking  is  accounted  for  by  seismic  coefficients  used  in 
the  lateral  force  formula,  which  correspond  to  the  soil 
profile  types  defined  in  the  latest  edition  of  the  Uniform 
Building  Code  (ICBO,  1997).  This  revision  also  contains 
a  "near-source"  factor  that  takes  into  consideration  ef- 
fects of  the  proximity  to  nearby  earthquake  source 
ruptures  on  shaking.  Maps  of  known  active  fault  near- 
-source  zones  have  been  prepared  for  use  with  the  1997 
UBC  (ICBO,  1998).  The  advisory  committee  believes 


that,  given  the  current  understanding  of  the  effects  of 
geologic  materials  and  structure  on  earthquake  ground 
motions,  there  would  be  no  benefit  in  establishing  "am- 
plified shaking  hazard  zones"  for  purposes  of  design  and 
construction.  The  purpose  of  the  Seismic  Hazard  Map- 
ping Act  is  to  identify  where  special  provisions,  beyond 
those  contained  in  the  UBC,  are  necessary  to  ensure 
public  safety.  This  need  has  not  been  recognized  for  the 
hazard  of  ground  shaking.  Design  provisions  contained 
in  the  UBC  are  believed  to  be  representative  of  current 
knowledge  and  capability  in  earthquake-resistant  de- 
sign. 
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Consideration  should  be  given  to  preparation  of  "infor- 
mational" maps  that  identify  where  soft-soil  profiles 
(type  St )  are  more  likely  to  be  found.  Similarly,  identi- 
fying areas  where  basin  structure  or  topography  may 
enhance  ground  shaking  or  where  an  aggregate  of  such 


adverse  conditions  within  near-source  zones  might  occur 
could  be  of  value  for  land-use  planning  purposes.  The 
development  and  utility  of  these  options  should  be  in- 
vestigated. 


LIQUEFACTION  HAZARD  ZONES 


INTRODUCTION 

California  Department  of  Conservation's  Division  of 
Mines  and  Geology  (DMG)  is  the  principal  state  agency 
charged  with  implementation  of  the  provisions  of  the 
1990  Seismic  Hazard  Mapping  Act.  These  recommen- 
dations are  developed  to  assist  DMG  in  mapping 
liquefaction  hazard  zones  (LQ-zones).  The  zones  estab- 
lish where  site-specific  geotechnical  investigations  must 
be  conducted  to  assess  liquefaction  potential  and,  if  re- 
quired, provide  a  technical  basis  to  mitigate  the 
liquefaction  hazard. 

Liquefaction  Hazard  Zoning  Criteria 

Liquefaction  hazard  zones  are  geographic  areas  meeting 
one  or  more  of  the  following  criteria: 

1.  Areas  known  to  have  experienced  liquefaction 
during  historic  earthquakes. 

Field  studies  following  past  earthquakes  indicate  lique- 
faction tends  to  recur  at  many  sites  during  successive 
earthquakes  (Youd,  1984).  There  are  many  published 
accounts  of  liquefaction  occurrences  and  the  areas  so 
delineated  should  be  included  in  the  LQ-zones. 

2.  All  areas  of  uncompacted  fills  containing  po- 
tentially liquefiable  material  that  are  saturated, 
nearly  saturated,  or  may  be  expected  to  become 
saturated. 

In  some  areas  there  has  been  a  practice  of  creating  us- 
able land  by  dumping  artificial  fill  on  tidal  flats  or  in  large 
deep  ravines.  Standard  geologic  criteria  are  of  little  use 
in  characterizing  soils  within  these  fills,  which  are  less  ho- 
mogeneous than  natural  deposits.  There  is  no  reason  to 
assume  lateral  stratification  in  these  fills  and  the  validity 


of  extrapolating  subsurface  data  is  questionable.  Evi- 
dence for  filling  can  be  found  using  maps  showing  old 
shorelines,  comparing  old  and  modern  topographic 
maps,  studying  logs  of  boreholes,  and  obtaining  reports 
or  original  plans  of  specific  projects  involving  reclaimed 
land.  These  areas  should  be  included  in  the  LQ-zones. 

3.  Areas  where  sufficient  existing  geotechnical 
data  and  analyses  indicate  that  the  soils  are  po- 
tentially liquefiable. 

The  vast  majority  of  liquefaction  hazard  areas  are  under- 
lain by  recently  deposited  sand  and/or  silty  sand.  These 
deposits  are  not  randomly  distributed,  but  occur  within  a 
narrow  range  of  sedimentary  and  hydrologic  environ- 
ments. Geologic  criteria  for  assessing  these  environ- 
ments are  commonly  used  to  delineate  bounds  of  suscep- 
tibility zones  evaluated  from  other  criteria,  such  as 
geotechnical  analysis  (Youd,  1991).  Ground  water  data 
should  be  compiled  from  well  logs  and  geotechnical 
boreholes.  Analysis  of  historical  aerial  photographs  may 
delineate  zones  of  flooding,  sediment  accumulation,  or 
evidence  of  historic  liquefaction.  The  Quaternary  geology 
should  be  mapped  and  age  estimates  assigned  based  on 
ages  reported  in  the  literature,  stratigraphic  relationships 
and  soil  profile  descriptions.  In  many  areas  of  Holocene 
and  Pleistocene  deposition,  geotechnical  and  hydrogeo- 
logic  data  are  available.  Geotechnical  investigation 
reports  with  Standard  Penetration  Test  (SPT)  and/or 
Cone  Penetration  Test  (CPT)  and  grain  size  distribution 
data  can  be  used  for  liquefaction  resistance  evaluations. 

For  sand  and  silty  sand,  there  are,  at  present,  two  accu- 
rate and  reliable  in-situ  approaches  available  for 
quantitative  evaluation  of  the  soil's  resistance  to  cyclic 
pore  pressure  generation  and/or  liquefaction.  These  are: 
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1)  correlations  and  analyses  based  on  in-situ  Standard 
Penetration  Test  (SPT)  D1586  (ASTM,  1990);  D6066- 
96el  (ASTM,  1999)  data,  and  2)  correlations  and 
analyses  based  on  in-situ  Cone  Penetration  Test  (CPT) 
ASTM  D3441  (ASTM,  1990)  data. 

Seed  and  others  (1984;  1985),  Seed  and  DeAlba 
(1986),  and  Youd  and  Idriss  (1997)  provide  guidelines 
for  performing  "standardized"  SPT.  They  also  provide 
correlations  for  converting  penetration  resistance  (ob- 
tained using  most  of  the  common  alternate  combina- 
tions of  equipment  and  procedures)  to  an  equivalent 
"standardized"  penetration  resistance  (Nj)^.  This  "stan- 
dardized" penetration  resistance  can  be  used  as  a  basis 
for  evaluating  liquefaction  resistance. 

Cone  penetration  test  (CPT)  tip  resistance  (qj  may  also 
be  used  as  a  basis  to  evaluate  liquefaction  resistance. 
This  is  done  either  by  empirical  comparison  between  q. 
data  and  case  histories  of  seismic  performance  (Olsen, 
1988)  or  by  converting  q^-values  to  "equivalent"  SPT  re- 
sistance and  use  of  correlations  between  (Nj)^  data  and 
case  histories  of  seismic  performance  (Robertson  and 
Campanella,  1985;  Seed  and  De  Alba,  1986;  Youd  and 
Idriss,  1997). 

In  addition  to  sandy  and  silty  soils,  some  gravelly  soils 
are  potentially  vulnerable  to  liquefaction.  At  present,  the 
best  available  technique  for  quantitative  evaluation  of 
the  liquefaction  resistance  of  this  type  of  deposit  involves 
correlation  and  analysis  based  on  in-situ  penetration  re- 
sistance measured  using  the  very  large  scale  Becker 
Hammer  system  (Harder,  1988;  1997). 

The  correlations  of  Seed  and  others  (1985),  as  updated 
in  Youd  and  Idriss  (1997),  and  the  (N,)^  data  can  be 
used  to  assess  liquefaction  susceptibility.  Because  geo- 
technical  analyses  are  usually  made  using  limited 
available  data,  the  susceptibly  zones  should  be  delineat- 
ed using  geologic  criteria.  Geologic  cross  sections,  tied 
to  boreholes  and/or  trenches,  should  be  constructed  for 
correlation  purposes.  The  units  characterized  by  geo- 
technical  analyses  should  be  correlated  with  surface  and 
subsurface  units  and  extrapolated  for  the  mapping 
project. 

Liquefaction  opportunity  is  a  measure,  expressed  in 
probabilistic  terms,  of  the  potential  for  ground  shaking 


strong  enough  to  generate  liquefaction.  Analyses  of  in- 
situ  liquefaction  resistance  require  assessment  of 
liquefaction  opportunity.  The  minimum  level  of  seismic 
excitation  to  delineate  LQ-zones  should  be  that  level  de- 
fined by  M7.5-weighted  peak  ground  surface 
acceleration  (PGA)  for  alluvial  soil  conditions  with  a 
10%  probability  of  exceedance  over  a  50-year  period. 

4.  Areas  where  existing  geotechnical  data  are  in- 
sufficient. 

In  areas  of  limited  or  no  geotechnical  data,  LQ-zones 
should  be  delineated  using  geologic  criteria  as  follows: 

(a)  Areas  containing  soil  deposits  of  late  Holocene  age 
(current  river  channels  and  their  historic  floodplains, 
marshes  and  estuaries),  where  the  M7.5-weighted  peak 
acceleration  that  has  a  10%  probability  of  being  exceed- 
ed in  50  years  is  greater  than  or  equal  to  0. 10  g  and  the 
historic  high  water  table  is  less  than  40  feet  below  the 
ground  surface;  or 

(b)  Areas  containing  soil  deposits  of  Holocene  age  (less 
than  11,000  years),  where  the  M7.5-weighted  peak  ac- 
celeration that  has  a  10%  probability  of  being  exceeded 
in  50  years  is  greater  than  or  equal  to  0.20  g  and  the 
historic  high  water  table  is  less  than  or  equal  to  30  feet 
below  the  ground  surface;  or 

(c)  Areas  containing  soil  deposits  of  latest  Pleistocene 
age  (between  11,000  and  15,000  years),  where  the 
M7.5-weighted  peak  acceleration  that  has  a  10%  proba- 
bility of  being  exceeded  in  50  years  is  greater  than  or 
equal  to  0.30  g  and  the  historic  high  water  table  is  less 
than  or  equal  to  20  feet  below  the  ground  surface. 

The  Quaternary  geology  may  be  taken  from  existing 
maps,  where  available,  and  hydrologic  data  should  be 
compiled. 

Candidate  Methods  for  Future  Development 

To  further  improve  delineation  of  liquefaction  zones  and 
strengthen  the  justification  for  geotechnical  site  investiga- 
tions, DMG  should  follow  the  development  of  methods 
based  on  quantifying  ground  deformation  associated 
with  the  occurrence  of  liquefaction.  Estimates  of  lique- 
faction potential  based  on  simplified  methods  are  known 
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to  be  conservative  with  regard  to  damage  potential.  Sur- 
face manifestation  of  liquefaction,  such  as  venting  of 
sand,  may  not  always  correlate  with  structural  damage, 
especially  when  only  a  small  fraction  of  the  soil  column 
liquefies  and  is  accompanied  by  little  or  no  settlement. 
Total  thickness  of  liquefiable  material  and  related  poten- 
tial for  significant  vertical  settlement  or  horizontal 


deformation  are  better  indicators  of  damage  potential. 
Improvements  in  generalized  measures  such  as  the  Liq- 
uefaction Potential  Index  (Iwasaki  et  al,  1982), 
Liquefaction  Severity  Index  (Youd  and  Perkins,  1987), 
and  displacement  from  lateral  spreading  (Bartlert  and 
Youd,  1995),  should  be  investigated  for  applicability  in 
delineating  seismic  hazard  zones  in  California. 


EARTHQUAKE-INDUCED  LANDSLIDE  HAZARD  ZONES 


INTRODUCTION 

The  California  Department  of  Conservation,  Division  of 
Mines  and  Geology  (DMG)  is  the  principal  state  agency 
charged  with  implementation  of  the  provisions  of  the 
1990  Seismic  Hazard  Mapping  Act.  These  recommenda- 
tions are  developed  to  assist  DMG  in  mapping 
earthquake-induced  landslide  hazard  zones. 

Landslide  Hazard  Zone  Criteria 

Earthquake-induced  landslide  hazard  zones  are  areas 
meeting  one  or  more  of  the  following  criteria: 

1.  Areas  known  to  have  experienced  earthquake- 
induced  slope  failure  during  historic  earthquakes. 

It  is  very  difficult,  if  not  impossible,  to  distinguish  earth- 
quake-induced slope  failures  from  landslides  triggered  by 
other  mechanisms  if  the  latest  movement  occurred  prior 
to  historic  observations.  Evidence  of  earthquake  trigger- 
ing for  large  pre-historic  landslides  tends  to  be  circum- 
stantial (for  example,  large  dormant  landslide  complexes 
are  often  near  active  faults),  and  the  shallow  disrupted 
landslides  (debris  or  soil  falls)  found  to  be  so  common  in 
historic  earthquakes  are  not  generally  preserved  in  the 
geologic  record.  However,  landslides  caused  by  some 
historic  earthquakes  in  California  have  been  well  docu- 
mented (Lawson,  1908;  Morton,  1975;  Harp  and  others, 
1984;  Spittler  and  Hart,  1990;  Harp  and  Jibson,  1995). 
Wherever  possible,  DMG  should  include  documented 
earthquake-triggered  landslides  within  zones  of  required 
investigation. 


2.  Areas  identified  as  having  past  landslide  move- 
ment, including  both  landslide  deposits  and 
source  areas. 

Steep  scarps  and  toe  areas  of  existing  landslides  often 
fail  in  moderate  to  large  earthquakes.  The  entire  mass 
of  existing  large  rotational  landslide  deposits  is  not  typi- 
cally reactivated  by  earthquake  shaking  (Keefer,  1984). 
However,  long-duration  earthquakes,  such  as  a  magni- 
tude 8+  earthquake  on  the  San  Andreas  Fault  in 
southern  or  northern  California,  could  reactivate  existing 
landslides  and  result  in  significant  damage  to  structures. 
Because  of  this  possibility,  existing  landslide  deposits 
and  their  source  areas  should  be  identified  and  included 
in  zones  of  required  investigation. 

An  inventory  of  all  landslides  should  be  prepared  for 
each  hazard  zone  map  area.  Al  existing  landslides,  in- 
cluding the  source  (scarp)  and  deposit,  should  be 
mapped  and  given  a  level  of  confidence  of  interpreta- 
tion. Landslides  identified  as  "definite"  or  "probable" 
should  be  added  to  the  geologic  strength  map  and 
should  always  be  included  in  zones  of  required  investiga- 
tion. Landslides  identified  as  "questionable,"  that  is, 
areas  having  geomorphic  features  that  are  probably  the 
result  of  other  causes  (e.g.,  stream  terraces)  and  would 
require  extensive  exploration  to  verify  a  landslide  origin, 
should  be  excluded  from  the  earthquake-induced  land- 
slide zones. 
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3.  Areas  where  DMG's  analyses  of  geologic  and 
geotechnical  data  indicate  that  the  geologic  mate- 
rials are  susceptible  to  earthquake-induced  slope 
failure. 

The  recommended  procedure  for  these  analyses  is  the 
Newmark  method  as  calibrated  by  McCrink  and  Real 
(1996),  described  below. 


Newmark  Method 

Currently,  the  most  advanced  method  for  mapping  re- 
gional earthquake-induced  landslide  hazards  is  based  on 
the  work  of  Newmark  (1965).  Newmark,  recognizing  the 
limitations  of  a  factor  of  safety  approach  to  dynamic 
slope  stability  analyses,  devised  a  method  of  estimating 
the  magnitude  of  ground  displacement  caused  by  a  given 
earthquake  ground  motion.  The  USGS  tested  New- 
mark's  method  on  a  landslide  triggered  by  the  1979 
Coyote  Lake  earthquake  (Wilson  and  Keefer,  1983),  and 
pioneered  the  application  of  the  Newmark  analysis  for 
mapping  earthquake-induced  landslide  hazard  potential 
in  San  Mateo  County  (Wieczorek  and  others,  1985). 

McCrink  and  Real  (1996)  calibrated  the  San  Mateo 
County  mapping  methodology  using  landslides  and  near- 
field  strong-motion  records  from  the  1989  Loma  Prieta 
earthquake.  They  also  developed  specific  procedures 
allowing  the  method  to  be  run  on  a  geographic  informa- 
tion system  (G1S).  Because  of  the  extensive  calibration 
and  validation  of  this  technique,  earthquake-induced 
landslide  hazard  zones  should  be  based  on  a  Newmark 
dynamic  displacement  analysis  using  the  parameters  and 
specific  approaches  that  have  been  developed  and  docu- 
mented by  McCrink  and  Real. 

The  following  paragraphs  briefly  describe  the  recom- 
mended analytical  procedure  developed  in  this 
calibration  study. 

Assumptions  in  the  Model 

In  order  to  efficiently  delineate  the  earthquake- induced 
landslide  zones  on  a  regional  basis,  the  following  as- 
sumptions and  simplifications  are  reasonable: 


•  The  failure  should  be  assumed  to  be  an  infi- 
nite-slope failure,  that  is,  a  relatively  shallow 
slide  that  has  a  failure  surface  parallel  to  the 
ground  surface. 

•  Only  unsaturated  slope  conditions  should  be 
considered. 

•  The  response  of  the  geologic  materials  to 
earthquake  shaking,  in  terms  of  landslide 
failure  potential,  should  be  characterized  by 
the  shear  strength  properties  of  the  geologic 
materials. 

Shear  Strength  Properties 

In  selecting  representative  shear  strength  properties  to 
characterize  geologic  materials,  DMG  should  use  the 
most  appropriate  combination  of  strength  parameters 
available  for  the  hazard  map  area.  The  calibration  study 
(McCrink  and  Real,  1996)  indicates  that  the  internal  an- 
gle of  friction  alone  is  adequate  for  regional  mapping  of 
earthquake-induced  slope  failure  potential.  Where 
appropriate,  DMG  should  identify  adverse  bedding  con- 
ditions (out-of-slope  bedding)  and  apply  shear  strength 
values  representing  the  weaker  materials  (such  as  shale 
interbeds  in  a  predominantly  sandstone  formation)  of  the 
mapped  geologic  unit.  If  geotechnical  shear  test  data  are 
sufficient  or  lacking  for  a  mapped  geologic  unit,  such  a 
unit  should  be  grouped  with  lithologically  and  strati- 
graphically  similar  units  for  which  shear  strength  data  are 
available.  Published  shear  strength  values  can  be  used  if 
necessary.  The  product  of  the  shear  strength  character- 
izations should  be  a  geologic  material  strength  map, 
wherein  the  areas  depicted  on  the  map  no  longer  repre- 
sent chronostratigraphic  "formations"  but  areas  of 
similar  shear  strength. 

Slope  Stability  Calculations  and  Factor  of 
Safety 

Slope  stability  calculations  using  the  infinite-slope  failure 
model  should  consist  of  first  calculating  a  static  factor  of 
safety,  followed  by  a  calculation  of  the  yield  acceleration 
from  Newmark' s  equation: 

a   =  (  FS  -  1  )g  sin  a 
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where  a  is  the  yield  acceleration  (the  horizontal  ground 
acceleration  required  to  cause  the  factor  of  safety  to 
equal  1 .0),  FS  is  the  factor  of  safety  from  the  static  sta- 
bility analysis,  g  is  the  acceleration  due  to  gravity,  and  is 
the  direction  of  movement  of  the  slide  mass,  in  degrees 
measured  from  the  horizontal,  when  displacement  is  ini- 
tiated (Newmark,  1965).  For  an  infinite-slope  failure 
model,  a  is  the  same  as  the  slope  angle. 

Earthquake  Ground  Motion 

Determination  of  anticipated  earthquake  shaking  for  the 
hazard  map  area  should  be  made  by  selecting  a  repre- 
sentative strong-motion  record  or  records,  based  on 
estimates  of  probabilistic  ground  motion  parameters  for 
levels  of  earthquake  shaking  having  a  10%  probability  of 
being  exceeded  in  50  years  (Petersen  and  others,  1996: 
Cramer  and  Petersen,  1996).  The  ground  motion  pa- 
rameters used  in  the  record  selection  should  include 
mode  magnitude,  mode  distance,  and  peak  acceleration. 

The  currently  recommended  procedure  calls  for  the  se- 
lected strong-motion  record  to  be  integrated  twice  for  a 
given  yield  acceleration  to  find  the  corresponding  New- 
mark  displacement.  This  process  should  be  repeated  for 
a  number  of  yield  accelerations  to  develop  a  mathemati- 
cal relationship  between  the  two  parameters.  The  yield 
acceleration  values  calculated  in  the  slope  stability  anal- 
yses should  be  correlated  with  Newmark  displacements 
estimated  from  the  strong- motion  record  to  prepare  a 
hazard  potential  map. 

Slope  Factors 

DMG  should  use  the  most  accurate  and  up-to-date  ter- 
rain data  available  to  derive  slope  and  aspect  maps. 
Digital  terrain  data  should  have  a  minimum  vertical  ac- 
curacy of  7  meters,  and  a  maximum  horizontal 
resolution  of  10  meters.  Acceptable  sources  of  terrain 
data  include  Level  2  digital  elevation  models  (DEMs) 
prepared  by  the  USGS,  terrain  data  derived  from  inter- 
ferometric  synthetic  aperture  radar,  photogrammetrically 
produced  terrain  data,  and  ground  survey  data.  The  se- 
lected terrain  data  sources  should  meet  or  exceed  the 
above  accuracy  and  resolution  requirements.  Slope  gra- 
dient and  slope  aspect  maps  prepared  from  the  digital 
terrain  data  should  be  generated  using  algorithms  most 
appropriate  for  the  terrain  data  used. 


Earthquake-Induced  Landslide  Potential 

An  earthquake-induced  landslide  potential  map  should 
be  prepared  by  combining  and  comparing  (overlay)  the 
geologic-material  strength  map  with  a  slope  gradient 
map.  Hazard  potential  criteria  for  the  hazard  maps 
should  be  based  on  the  amount  of  calculated  Newmark 
displacement  and  corresponding  slope  angle  for  each 
geologic  unit  caused  by  the  selected  strong-motion 
record:    "Very  Low"  would  correspond  to  displacements 
less  than  5  cm:  "Low"  potential  has  displacements  of  5 
cm  to  less  than  15  cm;  "Moderate"  potential  has  dis- 
placements of  15  cm  to  less  than  30  cm;  and  "High" 
potential  has  displacements  of  30  cm  or  greater.  On  the 
basis  of  the  calibration  study  (McCrink  and  Real,  1996), 
High,  Moderate  and  Low  levels  of  hazard  potential  (all 
areas  with  calculated  displacements  greater  than  5cm), 
should  be  included  within  the  landslide  zone  of  required 
investigation. 

Hazards  Not  Addressed 

Because  of  the  many  simplifying  assumptions  made 
when  applying  the  Newmark  analysis  to  regional  hazard 
mapping,  the  current  method  does  not  capture  all  types 
of  ground  failures  known  to  occur  during  earthquakes. 
Earthquake-generated  ground  failures  that  are  not  ad- 
dressed by  the  Newmark  method  include  those 
associated  with  ridge-top  spreading  and  shattered  ridges. 
Also,  run-out  areas  of  triggered  landslides  may  extend 
beyond  zone  boundaries  into  areas  outside  the  zone  of 
required  investigation.  The  potential  for  ground  failure 
resulting  from  liquefaction-induced  lateral  spreading  of 
alluvial  materials,  considered  by  some  to  be  a  form  of 
landsliding,  should  not  be  specifically  addressed  by  the 
earthquake-induced  landslide  hazard  zone  because  such 
hazards  are  to  be  included  in  the  LQ-zones. 

Candidate  Methods  for  Future  Development 

In  order  to  improve  the  accuracy  of  the  Newmark  meth- 
od in  capturing  all  appropriate  landslide-prone  areas, 
DMG  should  continue  to  refine  the  method.  From  re- 
cent earthquakes  we  know  that  ridge-top  spreading 
typically  occurs  along  strike-ridges,  and  that  shattered 
ridges  typically  occur  along  the  tops  of  high,  narrow  ridg- 
es.  DMG  should  use  this  knowledge  to  develop  models 
to  assess  the  potential  for  these  ground  failures  in  the  fu- 
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ture.  Methods  to  identify  rock  fall  and  debris  flow  run- 
out areas  should  also  be  investigated,  if  deemed 
adequate,  and  incorporated  into  future  zone  maps. 

In  addition  to  improving  the  current  Newmark  model, 
DMG,  in  cooperation  with  USGS,  should  continue  to  in- 
vestigate other  analytical  methods  that  might  be  useful 
in  zoning.  It  is  recommended  that  DMG  investigate  the 
applicability  of  two  analytical  methods  as  possible  alter- 
natives to  the  Newmark  method: 

1.  The  Multivariate  Method 

The  multivariate  method,  described  by  Carrara  and  oth- 
ers (1991),  uses  a  multivariate  statistical  procedure  in 
conjunction  with  GIS  techniques  to  model  landslide  haz- 
ards. In  this  method,  the  morphological,  geological  and 
vegetation  characteristics  for  slopes  are  analyzed  using  a 
stepwise  discriminant  analysis,  rating  the  characteristics 
in  terms  of  their  ability  to  discriminate  between  stable 
and  unstable  slopes.  The  method  does  not  specifically 
address  triggering  mechanisms  such  as  earthquakes  or 
rainfall,  but  holds  the  potential  to  identify  susceptible  ar- 
eas on  the  basis  of  past  performance  of  the  terrain  and 
other  characteristics. 


2.  The  Probabilistic  Slope  Stability  Method 

The  probabilistic  slope  stability  method  provides  a  sys- 
tematic and  quantitative  way  to  deal  with  uncertainties 
associated  with  soil  and  rock  spatial  variability,  geotech- 
nical  sampling  and  testing,  terrain  models,  and  earth- 
quake shaking.  Vanmarcke  (1976;  1980)  has  consid- 
ered the  basic  3-dimensional  stability  problem  in  a 
probabilistic  framework  for  man-made  embankments 
and  natural  slopes.  The  probabilistic  approach  has  the 
advantage  of  being  able  to  address  the  spatial  variability 
of  strength  parameters  and  ground-water  conditions, 
and  may  allow  for  the  easy  integration  of  probabilistic 
ground  motion  estimates.  The  USGS  is  evaluating  a 
form  of  probabilistic  earthquake  slope  stability  in  south- 
ern California  using  a  Newmark  displacement  model 
and  ground  motion  characterized  by  Arias  intensity  (Jib- 
son  and  others,  1998). 

The  multivariate  and  probabilistic  methods,  used  in  full 
or  in  part,  may  prove  suitable  as  possible  alternative  ap- 
proaches to  earthquake-induced  landslide  hazard 
mapping.  These  methods  are  not  currently  well  devel- 
oped for  regional  mapping  purposes,  and  calibration 
studies  will  need  to  be  conducted.  However,  some  or  all 
of  the  procedures  could  be  applied  to  more  accurately 
and  cost-effectively  delineate  earthquake- induced  land- 
slide hazard  zones. 
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